function [Q] = f3v2(lambda)

p0 = 2;
N = 100;
g = 9.8;

h = p0/N;

Q = 0;
intGamma = 0;
for i = 1:N
	intGamma = intGamma - (feval(@gamma, (i-1)*h)+feval(@gamma, i*h))*h/2;
	Q = Q + h/(sqrt(lambda + 2*intGamma));
end
Q = 2*g*Q + lambda;


end
